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Abstract 

Method of polarized semi-inclusive deep inelastic scattering (SIDIS) data analysis in 
the next to leading order (NLO) QCD is developed. Within the method one first directly 
extracts in NLO few first truncated (available to measurement) Mellin moments of the 
quark helicity distributions. Second, using these moments as an input to the proposed 
modification of the Jacobi polynomial expansion method (MJEM), one eventually recon- 
structs the local quark helicity distributions themselves. All numerical tests demonstrate 
that MJEM allows us to reproduce with the high precision the input local distributions 
even inside the narrow Bjorken x region accessible for experiment. It is of importance that 
only four first input moments are sufficient to achieve a good quality of reconstruction. 
The application of the method to the simulated SIDIS data on the pion production is 
considered. The obtained results encourage one that the proposed NLO method can be 
successfully applied to the SIDIS data analysis. The analysis of HERMES data on pion 
production is performed. To this end the pion difference asymmetries are constructed from 
the measured by HERMES standard semi-inclusive spin asymmetries. The LO results of 
the valence distribution reconstruction are in a good accordance with the respective lead- 
ing order SMC and HERMES results, while the NLO results are in agreement with the 
existing NLO parametrizations on these quantities. 

Introduction 

■ One of very important topics in the modern high energy physics is the investigation of the 
^ partonic spin structure of nucleon. In this connection, nowadays, there is a huge growth of 
^ ', interest to the SIDIS experiments with longitudinally polarized beam and target such as SMC 
- - ' P, HERMES 0, COMPASS ^. It is of importance that the SIDIS experiments, where 
one identifies the hadron in the final state, provide us with the additional information on 
the partonic spin structure in comparison with the usual DIS experiments. Namely, on the 
contrary to the DIS data, the SIDIS data allows us to extract the sea and valence quark 
helicity distributions in separation. 

At the same time it is argued (see, for example, Ref. [4 ) that to obtain the reliable 
distributions at relatively low average available to the modern SIDIS experiments^, the 
leading order (LO) analysis is not sufficient and NLO analysis is necessary. In Ref. |S] it 
was proposed the procedure allowing the direct extraction from the SIDIS data of the first 
moments (truncated to the accessible for measurement Bjorken x region) of the quark helicity 
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"^For example, HERMES data on semi-inclusive asymmetries is obtained at Q^noan = 2.5 GeV^ 
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distributions in NLO QCD. However, in spite of the special importance of the first moments^, 
it is certainly very desirable to have the procedure of reconstruction in NLO QCD of the 
polarized densities themselves. At the same time, it is extremely difficult to extract the local 
in Bjorken x (xb) distributions directly, because of the double convolution product entering 
the NLO QCD expressions for the semi-inclusive asymmetries^ (see Ref. j^] and references 
therein). Fortunately, operating just as in Ref. one can directly extract not only the first 
moments, but the Mellin moments of any required order. Using the truncated moments of 
parton distribution functions (PDFs) and applying the modified Jacobi polynomial expansion 
method (MJEM) proposed in Ref. [Hj one can reconstruct PDFs themselves in the entire 
accessible for measurement xb region. In the brief letter [Hj MJEM was tested using only the 
simple numerical (idealized) examples, where the exact values of the input moments entering 
MJEM are known (see Section 2). However, in the conditions of the experiment we have at our 
disposal only rather small number of the measured asymmetry values (one point for each bin 
with the rather wide bin widths at the middle and large xb)- Thus, extracting the moments 
from the measured asymmetries one calculates the integrals over x b using rather smalF number 
of points. So, because of this problem, even for the data obtained with very high precision (small 
errors) the extracted moments always suffer from the deviation from their true values. In this 
extended paper we investigate this problem in detail (see Section 3). In Section 3 MJEM is 
tested with the simulations corresponding to the kinematics of the HERMES experiment, where 
the accessible xb region is the most narrow in comparison with SMC and COMPASS regions. 

After the testing the proposed method is applied to the NLO QCD analysis of the HERMES 
data (Section 4). Notice that although the method is quite general, within this paper we deal 
only with the SIDIS data on the pion production. The point is that here we would like first of 
all to see how well the method itself works. That is why, for a moment, we do not like to deal 
with the such poorly known objects as and Dg fragmentation functions which additionally 
introduce the big uncertainties in the analysis results. For example, the analysis performed in 
Ref. [7j shows that the different choices of parametrizations for these fragmentation functions 
lead to the strong disagreement in the obtained results (about 30% for valence quarks and about 
100% for sea quarks). From this point of view the most attractive objects are the difference 
asymmetries jH] (for details on difference asymmetry in NLO see [01111 and references therein), 
where the fragmentation functions are cancel out in LO, while in NLO the difference asymmetry 
has only weak dependence of the difference of the favored and unfavored pion fragmentation 
functions (known with a good precision). In Section 4 the pion difference asymmetries are 
constructed from the measured by HERMES standard semi-inclusive spin asymmetries. Using 
as a starting point the constructed in such a way difference asymmetries, the (preliminary) 
reconstruction in NLO QCD of the valence PDFs from the HERMES data is performed. 

^Let us recall that namely these quantities, first moments, are of the most importance for solution of the 
proton spin puzzle because namely these quantities compose the nucleon spin 

^So, on the contrary to LO, where direct extraction of PDFs is possible, it seems at first sight that dealing 
with SIDIS asymmetries in NLO one can not avoid some fitting procedure. However, the modern world SIDIS 
data provide us by the rather small number of points for the measured asymmetries (and, besides, they suffer 
from the large statistical errors). Thus, purely semi-inclusive data very weakly constrains the large number of 
fit parameters entering NLO analysis (for example, twenty free parameters are used in Ref. 0). At the same 
time, the addition of DIS data in analysis can not help us to solve the main task of SIDIS - to extract the 
valence, sea and strange PDFs in separation. 

■^For example, HERMES used 9 bins for the region 0.023 < a; < 0.6 and COMPASS used 12 bins for the 
region 0.003 <x < 0.7. 
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2 MJEM and the usual JEM in comparison. Numerical 
tests. 



In this section we, for the sake of selfconsistence and clarity, represent in more detail the results 
ofRef. 0. 

There exist several methods allowing to reconstruct the local in xb quantities (like structure 
functions, polarized and unpolarized quark distributions, etc) knowing only finite number of 
numerical values of their Mellin moments. All of them use the expansion of the local quantity 
in the series over the orthogonal polynomials (Bernstein, Laguerre, Legendre, Jacobi) - see 
Ref. PU] and references therein. The most successful in applications (reconstruction of the 
local distributions from the evolved with GLAP moments and investigation of Aqcd) occurred 
the Jacobi polynomial expansion method (JEM) proposed in the pioneer work by Parisi and 
Sourlas [TT^ and elaborated^ in Refs. [12j and jlOj . 

The local in xb functions (structure functions or quark distributions) are expanded in the 
double series over the Jacobi polynomials and Mellin moments (see the Appendix): 

F{x) ^ F^_(x) = j-'^\x) yH:: 

^Y^U^tf^m + ii (1) 

where uj^°''^\x) = x^{l — x)"' and N^ax is the number of moments left^ in the expansion. For 
what follows it is of importance that the moments entering Eq. are the full moments, i.e., 
the integrals over the entire Bjorken x region < x < 1: 

M[j] = dxx^-^F{x). (2) 
Jo 

Till now nobody investigated the question of applicability of JEM to the rather narrow xb 
region available to the modern polarized SIDIS experiments. So, let us try to apply JEM to the 
reconstruction of Auy{x) and Adv{x) in the rather narrow xb region^" a = 0.023 < x < b = 0.6 
available to HERMES, and to investigate is it possible to safely replace the full moments (0) by 
the truncated moments. To this end we perform the simple test. We choose^^ GRSV2000NLO 
(symmetric sea) parametrization ^1] at = 2.5 GeV^. Integrating the parametrizations on 
Auv and Ady over the HERMES xb region we calculate twelve truncated moments given by 
(c.f. Eq. (Uni) below) 

M'[j] ^ M[;,] [j] ^ f dxx^-'Fix), (3) 

Ja 

where we put F{x) = Auy(x) or F{x) = Adv{x) and choose a = 0.023, 6 = 0.6 . Substituting 
these moments in the expansion Eq. (0) with N^ax = 12, we look for optimal values of 

^JEM with respect to polarized quark densities was first applied in Ref. J3| 

^Expansion |^ becomes exact when Nmax — * oo. However, the advantage of JEM is that even truncated 
series with the small number of used moments N„iax and properly fixed parameters a, (3 gives the good results 
(see, for example, [T7)] 1 

^°We choose here the most narrow HERMES xb region where the difference between JEM and its modification 
MJEM (see below) application becomes especially impressive. However, even with the more wide accessible Xb 
region (for example, COMPASS region [3]) it is of importance to avoid the additional systematical errors caused 
by the replacement of the full (unaccessible) moments in JEM Q by the accessible truncated moments. 
Certainly, one can choose for testing any other parametrization. 
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parameters a and /3 corresponding to the minimal deviation of reconstructed curves for Amv(x) 
and IS.dy[x) from the input (reference) curves corresponding to input parametrization. To find 
these optimal values Oopt and /3opt we use the program MINUIT [12]. The results are presented 
in Fig. n Looking at Fig. one can see that the curves strongly differ from each other even 
for the high number of used moments Nmax = 12. 

Thus, the substitution of truncated moments instead of exact ones in the expansion ^ 
is a rather crude approximation at least for HERMES xb region. Fortunately it is possible 
to modify the standard JEM in a such way that new series contains the truncated moments 
instead of the full ones. The new expansion looks as (see the Appendix) 



X — a\f^ f X — a 



h — a J \ b — a 



X E ©i"'^^ ( ^^ ) E cir^77-4uTT E luT^^m + l](-a)-', (4) 



where we use the notation Eq. Q for the moments truncated to accessible for measurement xb 
region. It is of great importance that now in the expansion enter not the full (unavailable) but 
the truncated (accessible) moments. Thus, having at our disposal few first truncated moments 
extracted in NLO QCD (Eqs. below), and applying MJEM, Eq. (jH), one can reconstruct 
the local distributions in the accessible for measurement xb region. 

To proceed let us clarify the important question about the boundary distortions. The 
deviations of reconstructed with MJEM, Eq. (0]), -F/v„„^ from F near the boundary points are 
unavoidable since MJEM is correctly defined in the entire region (a, b) except for the small 
vicinities of boundary points (see the Appendix). Fortunately, -F/v„„^ and F are in very good 
agreement in the practically entire accessible xb region, while the boundary distortions are 
easily identified and controlled since they are very sharp and hold in very small vicinities of 
the boundary points (see Figs. |2H below). In this section we, for clarity, explicitly show these 
distortions in all figures. In the next sections the all such distortions will be just cutted off. 

Let us check how well MJEM works. To this end let us repeat the simple exercises with 
reconstruction of the known GRSV2000NLO (symmetric sea) parametrization and compare 
the results of Auv{x) and Adv{x) reconstruction with the usual JEM and with the proposed 
MJEM. To control the quality of reconstruction we introduce the parameter^^ 

la ^-^l -^reconstructedi-^^ -^re/erence(''') I ]^QQ0/ 

eference \ 

where Freferencei.^) corrcspouds to the input parametrization and Freconstructed{.x) = F^^^^{x) 
in Eq. (^. We first perform the reconstruction with very high number of moments Nmax = 
12 and then with the small number Nmax = 4. Notice that the last choice Nmax = 4 is 
especially important because of peculiarities of the data on asymmetries provided by the SIDIS 
experiments. Indeed, the number of used moments should be as small as possible because first, 
the relative error |5(M'[j])/M'[j]| on M'[j] becomes higher with increase of j and second, the 
high moments become very sensitive to the replacement of integration by the sum over the 



^^Calculating i/ we just cut off the boundary distortions which hold for MJEM in the small vicinities of the 
boundary points (see the Appendix), and decrease the integration region, respectively. To be more precise, 
one can apply after cutting some extrapolation to the boundary points. However, the practice shows that the 
results on v calculation are practically insensitive to the way of extrapolation since the widths of the boundary 
distortion regions are very small (about 10^"^). 
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bins. The results of A-Uy(x) and Adv{x) reconstruction with MJEM at N^ax = 12 and with 
apphcation of both JEM and MJEM (in comparison) at Nmax = 4, are presented in Figs. El 
and El It is seen (see Fig. that for Nmax = 12 MJEM, on the contrary to the usual JEM 
(see Fig. Q), gives the excellent agreement between the reference and reconstructed curves. 
In the case Nmax = 4 the difference in quality of reconstruction between JEM and MJEM 
(see Fig. El) becomes especially impressive. While for standard JEM the reconstructed and 
reference curves strongly differ from each other, the respective curves for MJEM are in a good 
agreement. Thus, one can conclude that dealing with the truncated, available to measurement, 
xb region one should apply the proposed modified JEM to obtain the reliable results on the 
local distributions. 

Until now we looked for the optimal values of parameters a and (3 entering MJEM using 
explicit form of the reference curve (input parametrization) . Certainly, in reality we have no 
any reference curve to be used for optimization. However, one can extract from the data in 
NLO QCD the first few moments (Eqs. (jHl) below). Thus, we need some criterion of MJEM 
optimization which would use for optimization of a and (3 only the known (extracted) moments 
entering MJEM. 

On the first sight it seems to be natural to find the optimal values of a and (3 minimizing the 
difference of reconstructed with MJEM and input^^ {entering MJEM (01)) moments. However, 
it is easy to prove (see the Appendix) that this difference is equal to zero: 



reconstructed 



, n < Nmax, (6) 

input 



i.e. all reconstructed moments with n < Nmax are identically equal to the respective input 
moments for any a and (3. Fortunately, we can use for comparison the reference "twice- 
truncated" moments 

M"[n] = [n] = dx a;"-^F(x) {a < a + a' < b - b' < b), (7) 

' J a+a' 

i.e. the integrals over the region less than the integration region (a, b) for the "once-truncated" 
moments M^^^ entering MJEM The respective optimization criterion can be written in 
the form 



max 

E 

.,=0 



-^^(reconstructed) Ijl -^^(reference) Ijl 



mm. 



The "twice truncated" reference moments should be extracted in NLO QCD from the data in 
the same way as the input (entering MJEM (Q) "once truncated" moments. In reality one 
can obtain "twice-truncated" moments using Eqs. © (below) and removing, for example, first 
and/or last bin from the sum in Eq. ()14j) (below). 

Let us now check how well the optimization criterion ^ works. To this end we again perform 
the simple numerical test. We choose GRSV2000NLO parametrization at = 2.5 GeV^ with 
both broken and symmetric sea scenarios. We then calculate four first "once-truncated" and 
four first "twice-truncated" moments defined by Eqs. Q and and then substitute them in 
Eq. (jU and the optimization criterion (jSl), respectively. To find the optimal values of a and 

^•^In practice one should reconstruct these input moments from the data using Eqs. © (below). The reference 
"twice-truncated" moments (UJ should be reconstructed from the data in the same way. 
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j3 we use the MINUIT program. The results are presented in Fig. |3J It is seen that the 
optimization criterion works well for both symmetric and broken sea scenarios. 

Thus, the performed numerical tests show that the proposed modification of the Jacobi 
polynomial expansion method allows to reconstruct with a high precision the quark helicity 
distributions in the accessible for measurement xb region. 



3 Reconstruction of the valence quark helicity distribu- 
tions from the simulated data. 

In this section the proposed NLO QCD method will be applied to the simulated data. The 
simulations give us a good tool for testing of the method since here one knows in advance 
the answer to be found - the reference parametrization entering the generator as an input. 
At the same time the properly performed simulations (i.e., corresponding to the experimental 
statistics, binning and kinematical cuts) allow us finally adapt the method for application to 
the real experimental data. 

Let us first investigate the peculiarities of the nth moments extraction in the conditions of the 
real experimental binning (rather small number of bins covering the accessible for measurement 
Xb region). 

The simple extension of the procedure proposed in Ref. [S] gives for the n-th moments 
A„g = Jq^ dx x'^~^q{x) of the valence distributions the equations 



1 ^1") + A:'^ 



5L 



Andi 



{n)l 



-(n)2 



1 44"^ 



5L 



Ai! 



(n)l 



L 



(n)2 



(9) 



Here the notation is absolutely analogous to one used in Ref. 



4"' - £ 



n — 1 ATT^—TV 



dx 



{Auv — dy) 



X dzh[l + ®—Cqq®\{Di 



An) 
•^d — 



dx x^'-^Af-^' 



a. 



X / dZh[l + (^ — Cqq®\{Di 



A), 

{Uy + dy) 
D2), 



(10) 



where Di {D2) is favored (unfavored) pion fragmentation function, the quantities L(n)2 
are defined as 



L(n)l 


= ^(n 


L{n)2 


T TT^ 

= ^(n 


Th 


JZ 



^(n)u 
^{n)d 



^{n)d 

TTT- 



T-K 

^{n)d-> 

J 7r^ 
^(n)u^ 



LU, - 11 D^M) + i £ ^ ^nC,,{z')D^/-^ 



(12) 



where 



AnCqq{z) 



dx ^ACgg(X, z) 



are the nth moments of the polarized Wilson coefficients ACqq{x,z) entering the NLO ex- 
pressions for the difference asymmetries ^^^"'^ (the respective experimental expressions via 
counting rates are given by Eq. (fT9|) - see below): 

{AAuy - My) jl dZh[l + ®^ACqq®\{Di - D2) 

{Auv - dv) /i dz,,[l + ®^Cqq®\{D^ - D^) ' 



{Auy + Ady) !l dZh[l + ®^ACqq®\{D^ - D2] 
{Uy + dy) JzdZh[l + - D2) 



It should be noticed that in reality one can measure the asymmetries only in the restricted 
Bjorken x region a < a; < 6, so that the approximate equations for the truncated moments (c.f. 
Eq. m 



A'^q = M'[n] = f dxx^'-^Aqi^x) 

J a 



(13) 



of the valence distributions have the form Eq. with the replacement of the full integrals in 
Eq. (fTn|) by the sums over bins covering the accessible region a < x < b: 



i=l 



Uuy - dy){x,) / dZh[l + ®T^Cqq®\{D^ - D2), (14) 



and analogously for A^j^^ 

The approximation to Eq. (fT(H) given by Eq. p4|l is based on the assumption that all 
integrated quantities are the constants^^ within each bin. This is well-known "middle point" 
numerical integration method. 

However, it seems that there is a way to improve this approximation having in mind the real 
experimental situation. The point is that the reality of an experiment compel us to approximate 
by the constant within the bin only the measured quantity (difference asymmetry here), that 
can be written as 



A n'^ —71" 



£ Al-'-^'iixi)) e{x - x,.i)d{x, - x) 



1=1 



(15) 



where A'^'^ {{^i)) ^ is the mean value of asymmetry in ith bin, xq = a, x^^-^^ = b and 6{x) 
is the usual step function. At the same time, there is no any need to approximate by the 
constant another x-dependent quantities (unpolarized valence PDFs and Wilson coefficients 
here) entering the integrals over x^ as the known input. Thus, substituting Eq. ()15p in the 
initial integral equation Eq. (|TO|l we get (c.f. Eq. JH])) 



i=l 



Xi-i 



a.. 



rfxx"-^(4My - dv){x) / dZh[l + 0^Cqq®]{Di - D2). (16) 

Ztt 



Operating in such a way one puts the x-dependent measured quantity (difference asymmetries (a^i) 
here) to be equal its mean vahic in the ith bin, while the x-dependent rest is calculated in the point Xi = (xi). 
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Notice that the analogous way of the integral approximation was applied by the HERMES 
collaboration in Ref. [2^, where the moments A^g were reconstructed substituting extracted 
from the data ("measured") quantities (Ag/g)(xj) in the equation (see Eq. (46) in Ref. j2]) 



A'^q = dx 

"'0.023 j^]^ 



Aq 

— {{xi))9{x - Xi-i)9{xi - x) 

q 



x"-^g(x). (17) 



To see the advantage of Eq. (fTH|) application, let us compare it with the application of the 
integration procedure given by Eq. ()14|) . To this end we perform absolutely idealized LO test, 
where in each bin (we choose the HERMES binning) the value of asymmetry is directly calcu- 
lated from the given^^ parametrization on Auy and Ady using the theoretical LO expressions 
[Hj for the difference asymmetries 

^^n+-7r- _ 4Anv - Adv _ ^^+_^- Auy + Ady 



Auy — dy Uy + dy 

For simplicity, within this test we put {xi) = (xj — Xj_i)/2 {i = 1, . . . ,9; Xq = 0.023, Xg = 0.6), 
so that reconstructed with Eq. (|T8|l values of Auy{{xi)) and Ady{{xi)) exactly coincide with 
the respective input parametrization values in the points (xi) - see Fig El Now we calculate 
four first moments using reduced to LO Eqs. and the integration procedures given by Eqs. 
()14|) and ()16|) and then we apply^^ MJEM to both sets of the obtained moments. Looking at 
Fig. |31one can see that reconstructed in this way curves strongly differ from each other, and 
the curve obtained with application of the integration procedure given by Eq. (|T6|) is in much 
better agreement with the input (reference) parametrization. 

Thus, following the results of just performed test, from now on we will use namely Eq. (fTH|) 
performing the moments calculations. 

Let us now perform LO and NLO analysis of the simulated SIDIS data on 7r+ and vr^ 
production with both proton and deutron targets. To this end we use the PEPSI generator of 
polarized events fH]- The conditions of simulations are presented in Tableland correspond to 
the HERMES kinematics. Let us stress that all the cuts on Q^, xp, and in Tableware 
the standard physical^^ cuts applied by SMC, HERMES and COMPASS. The statistics 3 ■ 10^ 
in Tabled is the total number of DIS events for both proton and deutron targets and for both 
longitudinal polarizations. 



Table 1: Simulation conditions. Here xb and xp are the Bjorken and Feynman variables, respectively, Zh is 
the standard hadronic variable and W is the invariant mass of the final hadronic state. 



Elepton 


Xb 


Xp 


Zh 


27.5 GeV 


0.023 <xb < 0.6 


xp > 0.1 


Zh> Z = 0.2 




g2 


^ mean 


Events 


> 10 GeV^ 


Q2 > 1 GeV^ 


2.4 Gel/2 


3 ■ 10^ 



^^We choose for illustration GRSV2000LO (symmetric sea) parametrization. At the same time, it is easy to 
check that absolutely the same same picture holds for any other parametrization. 

^^Here we find aopt and Popt values requiring the minimal deviation of reconstructed with MJEM Auv{{xi)) 
and Adv{{xi)) from the reference parametrization values at the points (xi). 

^^For example, the important cut on invariant mass > 10 GeV^ is applied by these collaborations to 
exclude the events coming from the resonance region. 
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Using the simulated data we construct the difference asymmetries (see Ref. for details) 

1 (iVf/ - AT-- )Ltt - (iVf/ - Nf )Ln 



ATT ' - 



z PePtID (iVf + - iVff )LTt + (iVf + - iVff )Lu 



where A^fjQD are the counting rates integrated over in the region Z = Q.2 < < 1, 
-^Ti{TT) ~ ^Ti{TT)/^Ti{TT) luminosities, and the quantities p^, p^, / are equal to unity in 

the conditions of simulations with PEPSI. 

We first perform the LO analysis of the simulated difference asymmetries. The important 
peculiarity of LO analysis is that in this case one can perform the extraction of Auy and Ady 
in two ways. First is the direct extraction where one applies Eqs. fjl8|) in each bin - points 
with error bars in Fig. IHl The second method is the proposed one, where MJEM is applied to 
the LO extracted moments - dashed line in Fig. IHl The moments used in MJEM are extracted 
from the simulated difference asymmetries with application of reduced to LO Eqs. (p fTT^ . (fTH|) 
and are presented in Table El 



Table 2: Results for LO extracted truncated moments for the simulations with the entering PEPSI two 
different parametrizations: GRSV2000LO (symmetric sea) parametrization (top) and GRSV2000LO (broken 
sea) parametrization (bottom). For comparison the respective reference (obtained by direct integration of 
entering PEPSI input parametrizations) moments are also presented. 









A'Jv 




n 


Extracted 


Reference 


Extracted 


Reference 


1 


0.7042 ± 0.0124 


0.7176 


-0.2568 ± 0.0271 


-0.2618 


2 


0.1489 ± 0.0037 


0.1477 


-0.0439 ± 0.0079 


-0.0482 


3 


0.0467 ± 0.0016 


0.0457 


-0.0118 ± 0.0033 


-0.0135 


4 


0.0179 ± 0.0007 


0.0173 


-0.0041 ± 0.0015 


-0.0048 








Kdv 




n 


Extracted 


Reference 


Extracted 


Reference 


1 


0.5346 ± 0.0123 


0.5255 


-0.0952 ± 0.0274 


-0.1103 


2 


0.1318 ± 0.0036 


0.1282 


-0.0297 ± 0.0081 


-0.0331 


3 


0.0434 ± 0.0015 


0.0425 


-0.0098 ± 0.0034 


-0.0107 


4 


0.0167 ± 0.0007 


0.0166 


-0.0037 ± 0.0015 


-0.0039 



Looking at Fig. one can see that the input (reference) parametrization slightly deviates 
from both the directly extracted values of Awy and Ady and the reconstructed with MJEM 
curve. These deviations are unavoidable and are caused by the specific character of the events 
generation with PEPSI. Our experience shows that the asymmetries reconstructed from the 
generated events always slightly differ from the respective asymmetries calculated from the 
input parametrizations (with application of Eq. ()18|) in LO). One the other hand, comparing 
the directly extracted and the reconstructed with MJEM Auy and Ady, one can see that they 
are in a good agreement with each other. Thus, the performed LO testing encourage us that 
the proposed method of PDFs extraction could be successfully applied. 

Let us now clarify the important point concerning application of the optimization criterion 
Eq. (jS)) which we use to find the optimal values Uopt and (3opt of the entering MJEM parameters 
a and (3 (see Section 2 for details). All over the paper, applying the optimization criterion, we 
simultaneously use for each j (j = 1, . . . 4) in the sum two extracted twice-truncated moments. 
These moments correspond to two (sufficiently large and overlapping) integration regions, cov- 
ering respectively the bins from first to seven and from third to last ninth. We make such choice 
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because on the one hand, one should take into account in the criterion the whole accessible 
integration region, and, on the other hand, the "twice-truncated" moments should essentially 
differ from the "once truncated" moments for the well-working of the minimization procedure 
(see the respective discussion just after Eq. (jH))). 

Let us now perform the NLO analysis of the simulated data. We again use as an input two 
different parametrizations GRSV2000NLO (symmetric sea) and GRSV2000NLO (broken sea). 
The conditions of simulation are presented in Table H We first extract the truncated moments 
using Eqs. (P)- (fT^ and (fT^ . The results are presented in Table El Using these moments and 
applying MJEM, we reconstruct in NLO Auv{x) and Ady{x) with the results presented in Fig. 
13 Comparing Fig. [3 with Fig. IHl one can see that quality of reconstruction in NLO is not 
worse than the quality of LO reconstruction. The slight deviations of reconstructed and input 
curves as before (c.f. Fig. ^ are explained by the unavoidable deviations of the simulated with 
PEPSI asymmetries from their reference (corresponding to the input parametrization) values. 

Table 3: Results for NLO extracted truncated moments for the simulations with the entering PEPSI two 
different parametrizations: GRSV2000NLO (symmetric sea) parametrization (top) and GRSV2000NLO (broken 
sea) parametrization (bottom). For comparison the respective reference (obtained by direct integration of 
entering PEPSI input parametrizations) moments are also presented. 









Kdv 




n 


Extracted 


Reference 


Extracted 


Reference 


1 


0.7369 ± 0.0133 


0.7507 


-0.2577 ± 0.0293 


-0.2760 


2 


0.1507 ± 0.0039 


0.1545 


-0.0423 ± 0.0085 


-0.0490 


3 


0.0449 ± 0.0016 


0.0471 


-0.0109 ± 0.0033 


-0.0133 


4 


0.0163 ± 0.0007 


0.0176 


-0.0037 ± 0.0015 


-0.0045 








Kdv 




n 


Extracted 


Reference 


Extracted 


Reference 


1 


0.5860 ±0.0134 


0.5701 


-0.1045 ±0.0300 


-0.1137 


2 


0.1392 ±0.0039 


0.1381 


-0.0314 ±0.0088 


-0.0367 


3 


0.0433 ±0.0015 


0.0448 


-0.0101 ±0.0034 


-0.0121 


4 


0.0159 ±0.0007 


0.0172 


-0.0037 ±0.0016 


-0.0045 



The remark concerning very important peculiarity of application in NLO of the optimization 
criterion Eq. (jH} should be made here. The crucial point for the optimization criterion is the 
proper choice of the initial^^ values of a and (3. Indeed, the experience shows that if these 
initial values are too far away from the real aopt and f3opt to be found, then the MINUIT 
program can "fall" into some wrong local minimum and produce the false values of aopt and 
Popt- Fortunately in LO we can compare the reconstructed with MJEM curve with the reference 
(directly extracted) values of PDFs and unambiguously find the optimal values of aopt and Popt- 
On the other hand, it is natural to use aopt and jSopt obtained within LO analysis as the initial 
(starting) values for the application of optimization criterion Eq. (jH)) in NLO. The simulations 
demonstrate (see Fig. Ej) that with the such choice of initial a and (3 the minimization procedure 
performed in NLO analysis unambiguously finds the proper values of aopt and Popt- As a 
result the obtained with MJEM curves are in a good agreement with the input (reference) 
parametrizations. At the same time, looking at Figs. IHlandlHl one can see that the behavior 
of both LO and NLO extracted curves is in a good agreement with the respective behavior of 
the input (reference) parametrizations. 

^^These are starting values for the MIGRAD algorithm implemented in MINUIT package jl5| . 
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Thus, all performed in this section studies show that the proposed method can be success- 
fully applied to the polarized PDFs extraction in NLO QCD. 



4 NLO QCD analysis of the HERMES data on pion pro- 
duction. 

4.1 Construction of the difference asymmetries from the HERMES 
data on pion production 

Let us now apply the proposed method to the HERMES SIDIS data on the pion production. 
Within this paper we would like first of all to test the applicability of the method to the 
experimental data analysis. That is why, for a moment, we do not like to deal with the such 
poorly known objects as and fragmentation functions. As it was discussed before 

(see the Introduction), from this point of view the most attractive objects are the difference 
asymmetries. At the same time the difference asymmetries are still not constructed^^. So, let us 
apply a trick and express the difference asymmetry given by Eq. via the standard virtual 
photon SIDIS asymmetries 



/i7r± 



which were measured by HERMES 0. Namely, in each ith bin the difference asymmetries 
given by Eq. (|TI?|l can be rewritten as 

= f/^ - —j^ A--{xi), (20) 

where the quantity R^'' is defined as 

It is easy to see that the ratio Rt^~ can be rewritten as 

^+/_ ^ af;(x,) + <(x,) ^ aZ^poiix,) ^ 

and, thus, can be taken from the unpolarized SIDIS data. This relative quantity is well defined 
and extracted with the high precision object. We take its value from the LEPTO generator of 
unpolarized events ^Tj, which gives a good^° description of the fragmentation processes. The 
calculations show that the relative quantities R^^~ remarkably weakly depend on statistics of 



simulated events. Indeed, if one changes N^^^i from 10^ to 10® then only 1-3% deviation 

5+/- 



of i?j (in dependence of bin number) occurs. Nevertheless, to be more precise, extracting 
""^^At present the such analysis is performed by HERMES collaboration 

^^Dealing with LEPTO generator one should properly tune |2] the internal parameters of generator in order 
to achieve a proper description of the fragmentation process in the different experiments (HERMES here). 
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and Nl^^i 
condition 



491.000, iv-;„, 



the quantities Rj' for the HERMES statistics {N^otai = 117.000, N^^^^l = 82.000 for proton 



385.000 for deutron targets, respectively j2]) we preserve the 



total 



LEPTO 



total 



HERMES 



NLai = E,- , 



:22i 



performing the simulations with LEPTO. The results for K^^ are presented in Table 01 



Table 4: The obtained from the LEPTO generator results for the relative unpolarized quantity R^^ given 



by Eq. lETJ- 


proton target 


deutron target 






1 1.220 


1.150 


2 1.270 


1.201 


3 1.346 


1.229 


4 1.436 


1.274 


5 1.494 


1.315 


6 1.569 


1.350 


7 1.629 


1.407 


8 1.669 


1.444 


9 1.803 


1.556 



Thus, using in Eq. (p])) the results from Table El and HERMES results P] on A^^^ (see 
Tables XII and XIII in Ref. one easily constructs the difference asymmetries ^^^"'^ . The 
results are presented in Fig. 1101 

First, for the sake of testing (to check how well Eq. (j^IHl works), we reconstruct the valence 
PDFs in the leading order. In LO the equations for the difference asymmetries take the simple 
form given by Eq. (fTS|) . With these equations we reconstruct Auy and Ady using the results 
on difference asymmetries presented in Fig. ^1 The results are shown in Fig. ^2 , where we 
also plotted the respective results from Ref. ^ (obtained with the purity method). One can 
see that the results obtained with both procedures are in a good agreement. 

Let us now compare the LO extracted^^ moments (truncated to the HERMES xb region) 
obtained with application of Eq. (first line in the Table El) with the existing LO results of 
SMC and HERMES taken from the Table XI in Ref. [2]. One can see that the results obtained 
with application of Eq. ()20|) (i.e., from the difference asymmetries plotted in Fig. ^1) are in a 
good accordance with both HERMES and SMC results. 

Thus, the performed LO tests show that representation Eq. (f^n|) for the difference asym- 
metry can be successfully applied. 

Notice also that even the LO extraction of Auy and Ady from the difference asymmetries 
is interesting in itself as an alternative (complementary) possibility. Indeed, Eqs. ()18p are free 
from the rather badly known fragmentation functions and purities. 

^^As before (see Section 3), reconstructing the moments we apply the procedure of integration given by Eq. 
(|16|l since it gives more precise reconstruction of the local PDFs than the usually applied procedure given by 
Eq. (O. 



12 



Table 5: LO extracted truncated moments obtained from the difference asymmetries constructed with appH- 
cation of Eq. in comparison with the existing LO resuhs of SMC and HERMES coUaborations.The SMC 
moments are truncated to the HERMES xb region and are evolved to the HERMES Qmean = "2.5 GeV^ - see 
Table XI in Ref. [21. 





n 


1 


2 


3 


4 


This paper 


0.510 ± 0.110 


0.134±0.043 


0.048 ± 0.020 


0.020 ± 0.010 


HERMES 


0.603 ± 0.071 


0.144±0.014 


-/- 


-/- 


SMC 


0.614 ± 0.082 


0.152±0.016 


-/- 


-/- 




n 


1 


2 


3 


4 


This paper 


-0.280±0.146 


-0.074± 0.058 


-0.026 ± 0.026 


-0.011 ± 0.013 


HERMES 


-0.172±0.068 


-0.047± 0.012 


-/- 


-/- 


SMC 


-0.334±0.112 


-0.056± 0.026 


-/- 


-/- 



4.2 Reconstruction of the valence PDFs in NLO QCD 

Here, using the constructed difference asymmetries as a starting point and operating just as in 
Section 3, we will reconstruct in NLO QCD both the truncated Mellin moments of the valence 
PDFs and the local PDFs themselves. 

Let us first extract four first moments truncated to the HERMES xb region. The results are 
presented in Table IHl It is of importance that the proposed procedure allows us to extract the 
moments in NLO directly, without the commonly used assumptions like An = Ad = As = As 
(see, for example, Refs. O HE])- Notice that the first moments AiUy and Aidy are very 
important in themselves because namely the first moments compose the nucleon spin. At the 
same time all four moments presented in Table El are necessary for the reconstruction of the 
local PDFs with MJEM application. 

Table 6: NLO extracted truncated moments obtained from the difference asymmetries constructed with 
application of Eq. (|2()|l . 



n 


1 


2 


3 


4 


A>y 


0.555 ± 0.126 


0.134 ± 0.047 


0.047 ± 0.020 


0.019 ± 0.10 




-0.302 ± 0.173 


-0.076 ± 0.064 


-0.025 ± 0.027 


-0.010 ± 0.012 



Before application of MJEM in NLO, let us, for the sake of testing, reconstruct the local 
valence distributions in the leading order using LO moments from Table The results are 
presented in Fig. ^1 One can see that reconstructed with MJEM Eq. and optimization 
criterion Eq. (jH)) curve is in a good agreement with both HERMES results and with the results 
of direct LO extraction from the difference asymmetries constructed with application of Eq. 

dm). 

After the successful LO testing we apply MJEM in NLO QCD using the results for A'^uy 
and A'^dy from Table IHl The results are presented in Fig. ^1 where also, for comparison, the 
respective LO results are plotted. It is seen that the behavior of NLO and LO curves with 
respect to each other is in agreement with the predictions of existing parametrizations (see, for 
example, [Ej). 
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4.3 Corrections caused by dependence of asymmetries 

Until now we applied the approximation 

commonly used (see Refs. 0121) for analysis of the DIS and SIDIS asymmetries. This approx- 
imation is in a good agreement even with the COMPASS data (see Fig. 5 in Ref. [10]) and 
is especially suitable for the HERMES kinematics, where the "shoulder" in is rather small 
{IGeV^ < < lOGeV^; Q^an = '2.5 GeV^ - see Tables XII and XIII in Ref. PJ in compari- 
son with the SMC and COMPASS kinematics. Nevertheless, even for the HERMES kinematics 
we deal with, for more comprehensive analysis, it is useful to estimate the corrections caused 
by the weak dependence of the difference asymmetries. So, let us estimate the shifts in all 
four NLO moments caused by the respective shifts 

^i^p,d ^ ~ ^p,d ^ i^iiQmean) ~ ^p,d {^iiQl) (23) 

in the difference asymmetries. The most simple way to estimate 5iAp^^~'^ is to use the maximal 
number of the latest available NLO parametrizations. Namely, we approximate r.h.s of Eq. 

by the respective difference of "theoretical" asymmetries calculated with substitution of 
the different parametrizations to NLO equations for the difference asymmetries - Eqs. (14), 
(15) in Ref. 0. 

Adding the calculated in this way 6iAp^^~'^ to the initial experimental asymmetries A^^^'"" (xj, 

one estimates the evolved from to Q^ean asymmetries Ap^^'"^ {xi, Qmean)\Evoived- Using the 
obtained in such a way evolved asymmetries we extract the respective corrected moments of 
the valence PDFs A'^qvlcorrected repeating the procedure from Section 3. Then we compare 
the corrected moments A'^qvlcorrected with the respective moments A^gy from the previous 
section (obtained without corrections due to evolution) and calculate the respective shifts 
S{A'^qv) = 'AnQvlcorrected — ^'nlv ■ The results are presented in the Tables and |H1 where 
also the relative quantities (5(A^gv')/A^gy are presented. 

Notice that the considered procedure of the asymmetry evolution is quite similar to the 
procedure used by SMC for the rip(d) reconstruction (see Section V in Ref. pUj). 

It is of importance that we use for estimations the set of essentially different^^ NLO 
parametrizations and some of them, for example GRSV2000 (broken sea) and GRSV2000 (sym- 
metric sea), differ from each other very strongly. However, one can see from the Tables [3 |H1 that 
independently of the chosen parametrization the corrections for moments caused by evolution 
are very small (negligible in comparison with the statistical errors). To be precise, one can just 
include 5 {/S.[^qv)\ Average (see Tables m and IHI) in the systematical error. 

Let us now reconstruct the local valence PDFs applying MJEM to the corrected moments 
^'nlvlcorrected ffom the Tables [3 and |H1 The results are presented in Fig. HH One can see that 
the curves corresponding to the different ways of correction (different used parametrizations) 
are very close to each other. The averaged over the different corrections curve (dashed line) 
also very insignificantly differs from the initial (obtained without corrections) curve (solid line). 

Thus, the performed in this Section analysis demonstrates that at least for the HERMES 
kinematics we deal with here, the results are very insensitive to the corrections on the difference 
asymmetries caused by evolution. 

^^They correspond to the different sea scenarios (symmetric sea, weakly and strongly broken light quark sea), 
different details of calculations and different choice of ansatz for PDFs. 
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Table 7: NLO results for AJ^uy corrected due to evolution (top) together with the respective absolute and 
relative deviations from the uncorrected moments (bottom) . The corrections are estimated using seven different 
NLO parametrizations. The roman numbers I and II correspond to GRSV2000NLO parametrization for broken 
and symmetric sea scenarios, respectively. The roman numbers III and IV correspond to Ref. |21| for sets 
ii+ and ii- (symmetric and weakly broken sea scenarios). The rest of numbers V-VII corresponds to the NLO 
parametrizations from Ref. [18. (in the order of citation). 



^'n'^v\corrected 



n 


I 


II 


III 


IV 


V 


VI 


VII 


1 


0.5495 


0.5464 


0.5555 


0.5551 


0.5473 


0.5588 


0.5457 


2 


0.1364 


0.1367 


0.1378 


0.1377 


0.1368 


0.1387 


0.1363 


3 


0.0459 


0.0460 


0.0463 


0.0463 


0.0460 


0.0467 


0.0459 


4 


0.0182 


0.0182 


0.0183 


0.0183 


0.0182 


0.0185 


0.0182 








Average 


^fi^y 1 Corrected 







n 



1 0.5437 ± 0.1266 

2 0.1348 ± 0.0475 

3 0.0453 ± 0.0199 

4 0.0179 ± 0.0089 



n 


I 


II 


III 


IV 


V 


VI 


VII 


1 


-0.0054 


-0.0085 


0.0006 


0.0002 


-0.0077 


0.0039 


-0.0092 


2 


-0.0034 


-0.0031 


-0.0019 


-0.0021 


-0.0030 


-0.0010 


-0.0035 


3 


-0.0013 


-0.0012 


-0.0009 


-0.0009 


-0.0012 


-0.0005 


-0.0013 


4 


-0.0005 


-0.0005 


-0.0004 


-0.0004 


-0.0005 


-0.0002 


-0.0006 



^(A>y)/A>v^ (%) 



n 


I 


II 


III 


IV 


V 


VI 


VII 


1 


-0.97 


-1.53 


0.11 


0.04 


-1.38 


0.71 


-1.65 


2 


-2.42 


-2.21 


-1.37 


-1.47 


-2.12 


-0.73 


-2.47 


3 


-2.67 


-2.54 


-1.86 


-1.95 


-2.56 


-1.14 


-2.84 


4 


-2.62 


-2.62 


-2.08 


-2.14 


-2.72 


-1.28 


-2.94 



Average 5{/S.'^uv) 

n 

1 -0.0037 

2 -0.0026 

3 -0.0011 

4 -0.0004 



5 Conclusions and prospects 

Thus, in this paper the method of polarized SIDIS data analysis in NLO QCD is developed. 
The main peculiarity of the method is that its application is based on two subsequently applied 
procedures. First one directly extracts in NLO few first truncated (available to measurement) 
Mellin moments of the quark helicity distributions. The obtained at this stage results are very 
important and interesting in themselves. Indeed, the first moments are the main objects for 
understanding of the nucleon spin structure since they compose the nucleon spin. Second, using 
the obtained at first stage truncated moments as an input to the modification of the Jacobi 
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Table 8: NLO results for A'^dy corrected due to evolution (top) together with the respective absolute and 
relative deviations from the uncorrected moments (bottom). The corrections are estimated using seven different 
NLO parametrizations. The roman numbers I and II correspond to GRSV2000NLO parametrization for broken 
and symmetric sea scenarios, respectively. The roman numbers III and IV correspond to Ref. |21| for sets 
ii+ and ii- (symmetric and weakly broken sea scenarios). The rest of numbers V-VII corresponds to the NLO 
parametrizations from Ref. (18; (in the order of citation). 



A' ^ 1 

'-^n V 1 Corrected 


n 


I 


II 


III IV 


V 


VI 


VII 


1 


-0.3130 


-0.3091 


-0.3197 -0.3251 


-0.3062 


-0.3096 


-0.3064 


2 


-0.0778 


-0.0779 


-0.0811 -0.0822 


-0.0771 


-0.0780 


-0.0774 


3 


-0.0256 


-0.0256 


-0.0267 -0.0269 


-0.0254 


-0.0256 


-0.0255 


4 


-0.0096 


-0.0098 


-0.0102 -0.0102 


-0.0097 


-0.0097 


-0.0097 


Average A'^dylcorrected 



n 

1 -0.3127 ±0.1731 

2 -0.0788 ±0.0643 

3 -0.0259 ±0.0269 

4 -0.0099 ±0.0119 



n 


I 


II 


III 


IV 


V 


VI 


VII 


1 


-0.0114 


-0.0075 


-0.0181 


-0.0235 


-0.0046 


-0.0080 


-0.0048 


2 


-0.0015 


-0.0017 


-0.0048 


-0.0059 


-0.0008 


-0.0018 


-0.0012 


3 


-0.0003 


-0.0004 


-0.0014 


-0.0017 


-0.0001 


-0.0003 


-0.0002 


4 


-0.0001 


-0.0001 


-0.0005 


-0.0006 


-0.0000 


-0.0000 


-0.0000 








6iKdv)/Kdv 


(%) 






n 


I 


II 


III 


IV 


V 


VI 


VII 


1 


3.77 


2.48 


5.99 


7.78 


1.53 


2.65 


1.58 


2 


1.97 


2.20 


6.33 


7.72 


1.10 


2.33 


1.53 


3 


1.31 


1.47 


5.71 


6.74 


0.52 


1.35 


0.95 


4 


0.83 


0.83 


5.07 


5.89 


0.10 


0.52 


0.52 



Average 6{A'^d\ 



n 

1 -0.0111 

2 -0.0025 

3 -0.0007 

4 -0.0002 



polynomial expansion method, one eventually reconstructs the local quark helicity distributions 
in the accessible for measurement xb region. 

After successful testing we apply the proposed NLO method to the HERMES data on the 
pion production. To this end the pion difference asymmetries are constructed for both proton 
and deutron targets. To construct the difference asymmetries we use the HERMES data on 
the usual virtual photon SIDIS asymmetries, and, also, the well known quantity - ratio of 
unpolarized cross-sections for tt'^ and tt^ production which we take from the LEPTO generator 
of unpolarized events. With the constructed in such a way difference asymmetries the LO 
results of the valence distribution reconstruction are in a good accordance with the respective 
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leading order HERMES and SMC results, while the NLO results are in agreement with the 
existing NLO parametrizations on these quantities. Nevertheless, the obtained results should 
be considered as the rather preliminary since we construct the difference asymmetries in indirect 
way. At present the difference asymmetries are constructed by HERMES and are expected to 
be available in the nearest future. Certainly, it is very desirable to perform the NLO analysis 
of these directly constructed asymmetries and compare the results with the respective results 
presented in this paper. 

In this paper we apply the proposed method to the difference asymmetries only. Let us 
recall once again that essential advantage of these asymmetries in comparison with any other 
ones is the absence of fragmentation functions in LO and the weak dependence of well known 
difference of favored and unfavored pion fragmentation functions in NLO. So, since within this 
paper we mainly would like to investigate how well the method itself works, we, for a moment, 
prefer to deal namely with these very clean (from theoretical point of view) objects. In this 
connection it is of importance that the measurement of the difference asymmetries is one of the 
main topics of the physical program of E04-113 experiment planned at Jefferson Lab [22] ■ It 
is of importance that in this experiment the expected average is also rather small (about 
1GeV'^\ So, the NLO analysis in this experiment is also strongly required. 

Certainly, in the nearest future we will apply the method to NLO analysis of all measured 
in the SIDIS experiments asymmetries. In particular, it can allow us to extract in NLO so 
important quantity as the polarized strangeness in nucleon. Regretfully, the only existing 
today data on kaon production (HERMES experiment) suffers of large statistical errors and, 
besides, the accessible for measurement HERMES xb region is rather narrow. So, to obtain the 
reliable results on the such tiny quantity as Ais (as well as on Aim and Aid) it is necessary to 
perform the combined analysis, i.e., to analyze in NLO the combined data of SMC, HERMES, 
COMPASS and the planned E04-113 (Jefferson Lab) experiments. This is one of the main 
subjects of our future investigations. 
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Appendix 

JEM is the expansion of x-dependent function (structure function or quark density) in the 
series over Jacobi polynomials G^"'^''(x) orthogonal with weight a;*-"'^^(a;) = x^{\ — x)" (see 
-HO] for details): 

where 

M[j] = dxx^-^F{x) (A.2) 

JO 
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and 



dxJ'''''\x)el^'''\x)e^y\x) = 6„ 

The details on the Jacobi polynomials 



(A.3) 



et'\x) 



E 

j=0 



(A.4) 



can be found in Refs. JT] and jT^j. In practice one truncates the series ()A.1|) living in the 
expansion only finite number of moments Nmax ^ see Eq. (jij. The experience shows that 
even small Nmax gives good results. 

The idea of modified expansion is to reexpand F{x) in the series over the truncated moments 
^[a6]b] given by Eq. Q, performing the rescaling x ^ a + [b — a)x which compress the entire 
region [0, 1] to the truncated region [a, h]. To this end let us apply the following ansatz 



Fix) 



b — a 



b — a 



a oo 



n=0 



X — a 
b — a 



(A.5) 



and try to find the coefficients /„. Multiplying both parts of Eq. ()A.5|) by 6["'^''((a; — a)/(6 — a)), 
integrating over x in the limits [a, 6] and performing the replacement t = {x — a) / {b — a), one 
gets 



dxF{x)Q 



(",/3) 



X — a 
b — a 



n=0 



ih-a)Y.~fn r dtt^{l-tTQ^:'^\t)Qt'\t). (A.6) 



so that with the orthogonality condition Eq. ()A.3j) one obtains 



/„ = (6-a)-i / dxF{x)Q^: 



X — a 



b — a 



(A.7) 



Substituting Eq. ()A.7jl in the expansion ()A.5jl . and using Eq. ()A.4jl one eventually gets 



F{x) 



/x 


-ay 




-a. 


u 


— a J 




-a) 



E 



\k-l 



(A.8) 



where M'^j,][j] is given by Eq. 0. Truncating in the exact Eq. ()A.8|1 the infinite sum over 



n 



to the sum I]^=o^ one gets the approximate equation (jl)). 

Let us prove the important property^^ of the truncated moments reconstructed with MJEM. 
For any n < Nmax 



reconstructed 



input 



where 





= rrfxx"-^F(x) 




, J a 

mput 



reconstructed 



(A.9) 



(A. 10) 
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The proof of analogous property for the usual JEM can be found in ^HI 
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^Nmaxi.^) is the function reconstructed with apphcation of MJEM (jH), and N^ax + 1 is the 
number of the highest of moments used in Eq. 

To prove this statement we will need the inverse to ()A.4|) expansion 



(a,/3)o(«,/3), 



X] 



fc=0 



with the obvious property of d^nk^^ coefficients 



k=j 



(a,/3) 
nk kj 



Let US integrate Eq. Q over x in the limits [a, b] with weight {{x — a)/{b — a))". 



/ X ~ a 

ax 

a V — a 



LU 



(«,/3) 



, fx — a 

ax 

a \b — a 

^ max 



X 



(fe-a^E 



X — a 
b — a 



/ y m 



Ftv (x) 
X — a 



m=0 



j=0 



dzF{z) 



b — a 
z — a^^ 



b — a 



Using expansion ()A.11|) and orthogonality condition ()A.3|) . one easily gets 



(A.ll) 



(A.12) 



(A.13) 



( X — a 
dx 

a \b — a 



ax fT' 

E E ^km 

m=0 k=0 



(a,/3) (a,f3) 



(^nk E 



mj 



j=0 



dx Fix ) 



X — a\^ 



b — a 



(A. 14) 



It is obvious that at n < N^ax the Kronecker symbol 5km reduces the sum Ylm='o ^^e sum 
J2m=o- Thus, summing over m with 6km, using the identity 



n k 



n n 

E E — E E /j'^ 

k=Oj=0 j=Ok=j 



and applying Eq. ()A.12j) . one get eventually: 



b — a 



dx(-^ — -] F{x), n<N„ 
b — a J 



Setting n = in Eq. ()A.16j) one obtains 



^ fdxFN_S^)= f dxF{x)^M[^^,^[l\ 

J a J a 



reconstructed 

Putting then n = 1 in ()A.16jl and using ()A.17jl one gets 



input 



reconstructed 



input 



(A.15) 



(A.16) 



(A.17) 



(A.18) 
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Operating in this way for all n < N^ax, one arrives at the equality ()A.9|) to be proved. 

In conclusion, very important remark should be made here. Notice that ansatz ()A.5|) (as well 
as the expansion Eq. (0)) itself) is correctly defined inside the entire region (a, b) except for the 
small vicinities of boundary points (absolutely the same situation holds for the usual JEM, Eq. 
(jA.lll . applied to the quark distributions in the region (0, 1)). Thus, near the boundary points 
the deviations of reconstructed with MJEM function from its true values are unavoidable. 
Fortunately, all numerical examples (see section 2) show that input and reconstructed with 
MJEM functions are in very good agreement in the practically entire considered xb region, 
while the boundary distortions are easily identified and controlled since they are very sharp 
and hold in very small vicinities of the boundary points (see Figs. I21I1I). Thus, to reconstruct the 
curve near the boundary points one should just cut off these distortions and then extrapolate 
the rest to the boundaries of the considered xb region. 

To be precise, it should be also noticed that the all proofs given in the Appendix are rather 
formal because of the boundary distortions problem. All equations in the Appendix (like, for 
example, Eq. ()A.9|l ) become exact only when the distortions are cutted off and extrapolation 
to the boundaries is made. Fortunately, the practice show that the distortion regions are so 
small that the numerical results on the integrals over the entire region [a, b] are practically 
insensitive to the way of extrapolation, so that all equations in the Appendix are valid with a 
high numerical precision. 
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Figure 1: Results of Auv{x) [aopt = 8.189221, f3opt = -0.99000) and Adv{x) {aopt = -0.99000, Popt = 
—0.387196) reconstruction with the usual JEM. Solid line corresponds to the input (reference) parametrization 
GRSV2000NLO (symmetric sea). Dotted line corresponds to the distribution reconstructed with JEM for 
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Figure 2: Resuhs of Auy(x) {aopt = -0.827885, Popt = -0.011505) and Ady(x) {aopt = -0.989752, 
Popt — —0.012393) reconstruction with MJEM. Solid line corresponds to the input (reference) parametriza- 
tion GRSV2000NLO (symmetric sea). Dotted line corresponds to the distribution reconstructed with MJEM 

for Nmax = 12. 
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Figure 3: Results of the valence PDFs reconstruction with JEM and MJEM in comparison. Top part cor- 
responds to Auvix) [uopt = -0.99, Popt = 0.054010) and Mv{x) {uopt = 0.174096, /Sopt = 0.162567) recon- 
structed with the usual JEM. Bottom part corresponds to Auv{x) {aopt — —0.0025869, Popt — —0.071591) 
and Adv{x) {aopt = 0.110331, /3opt — —0.049255) reconstructed with MJEM. Solid lines correspond to in- 
put (reference) paramctrization GRSV2000NLO (symmetric sea). Dotted lines correspond to the distributions 
reconstructed with JEM (top) and MJEM (bottom). 




Figure 4: Results of Auy and Adv reconstruction for GRSV2000NLO paramctrization for both symmetric 
(top) and broken sea (bottom) scenarios. Solid line corresponds to the reference curve (input paramctriza- 
tion). Dotted line is reconstructed with MJEM and criterion ((HJ inside the accessible for measurement region 
([0.023,0.6] here). Optimal values of parameters for symmetric sea scenario for Auv are aopt = —0.15555, 
Popt — —0.097951 and for Ady are aopt = —0.002750, f3opt = —0.07190. Optimal values of parameters 
for broken sea scenario for Auy are aopt = —0.209346, Popt — 0.153417 and for Ady are aopt = 0.702699, 
= -0.293231. 
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Figure 5: (color online) Idealized LO testing of the integration procedures given by Eq. (|14|l and Eg. H16|l . The 
input parametrizations GRSV2000LO with symmetric sea scenario (top) and broken sea scenario (bottom) are 
used for direct calculation of asymmetries with Eq. I|18() . Solid line corresponds to the input parametrization. 
Closed circles correspond to the values of input parametrization in the middle of each bin. Broken line shows 
the way of integral approximation corresponding to application of Eq. (|14|l . Dashed line is obtained with 
MJEM and application of Eq. Hlt)|) for the moment calculation. Dot-dashed line is obtained with MJEM and 
application of Eq. H14|) for the moment calculation. Parameters vi and 1^2 are given by Eq. IS} and show the 
quality of reconstruction for the dashed and dot-dashed lines, respectively. 
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Figure 6: Results of the simulated difference asymmetry analysis in LO. GRSV2000LO parametrizations for 
symmetric sea (top) and broken sea (bottom) scenarios are used for simulations. Solid line corresponds to the 
input parametrization. Dashed line corresponds to the reconstructed with MJEM curve. Points with error bars 
correspond to direct extraction of the valence distributions with Eq. 1)18(1 . 
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Figure 7: Results of the simulated difference asymmetry analysis in NLO. GRSV2000NLO parametrizations 
for symmetric sea (top) and broken sea (bottom) scenarios are used for simulations. Solid line corresponds to 
the input parametrization. Dashed line corresponds to the reconstructed with MJEM curve. 
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Figure 8: Combined results of LO and NLO analysis (top) of the simulated difference asymmetries in com- 
parison with the respective LO and NLO versions of GRSV2000 (symmetric sea) parametrization (bottom). 
Dashed and solid lines correspond to LO and NLO curves, respectively. 
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Figure 9: Combined results of LO and NLO analysis (top) of the simulated difference asymmetries in compar- 
ison with the respective LO and NLO versions of GRSV2000 (broken sea) parametrization (bottom). Dashed 
and solid lines correspond to LO and NLO curves, respectively. 
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Figure 10: Difference asymmetries constructed with application of Eq. (|20|) . 
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Figure 1 1 : (color online) LO extraction of the valence PDFs from the difference asymmetries constructed with 
Eq. H2Q|I (up-oriented triangles) in comparison with the respective published HERMES results (down-oriented 
triangles). The HERMES results are shifted to the right for better visibility. 
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Figure 12: (color online) Different LO procedures of the valence PDFs extraction in comparison. Solid line 
and up-oriented triangles correspond to LO extraction with the proposed method and direct extraction with Eq. 
II18|I . respectively. The difference asymmetries constructed with application of Eq. H20() are used. Down-oriented 
triangles correspond to LO results of HERMES obtained with application of the purity method to the measured 
by HERMES usual virtual photon spin asymmetries. 
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Figure 13: Results of both LO and NLO analysis of the difference asymmetries constructed with Eq. H20|l . 
Solid and dashed lines correspond to NLO and LO results, respectively. 




Figure 14: (color online) NLO results obtained with application of MJEM to the moments extracted from 
the difference asymmetries corrected due to evolution. Dotted lines correspond to corrections estimated with 
the different parametrizations. Dashed line corresponds to the curve averaged over corrections. Solid line 
corresponds to reconstruction without corrections. 
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